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The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < pr < 
5.0 GeV/c) of electrons as a function of centrality in Au+Au collisions at ^/sjvjv = 200 GeV. 
Contributions from photon conversions and from light hadron decays, mainly Dalitz decays of tt° and 
n mesons, were removed. The resulting non-photonic electron spectra are primarily due to the semi- 
leptonic decays of hadrons carrying heavy quarks. Nuclear modification factors were determined by 
comparison to non-photonic electrons in p+p collisions. A significant suppression of electrons at 
high pr is observed in central Au+Au collisions, indicating substantial energy loss of heavy quarks. 

PACS numbers: 25.75.Dw 



It is well established that neutral pions and charged hadrons are strongly suppressed at high transverse mo- 
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mentum (pr) in high energy Au+Au collisions 0, 0, 0, 
0, Q . The suppression, which is absent in d+ Au colli- 
sions, implies that hard scattered partons traversing the 
medium created in Au+Au collisions experience consid- 
erable energy loss. Although high px suppression is ex- 
pected for charm quarks as well, their interaction with 
the medium has been predicted to be smaller than for 
light quarks, i.e. they should lose a lower fraction of 
their energy, as their large mass decreases the phase space 
available for gluon radiation, which is known as the "dead 
cone" effect p. If the medium is indeed less opaque to 
charm quarks they will also participate less in the collec- 
tive expansion of the medium, leading to a smaller elliptic 
flow strength vi [jj for particles carrying charm quarks 
compared to those solely composed of light quarks. Such 
medium effects should be even less pronounced for bot- 
tom than for charm quarks. 

The interaction of heavy quarks with the medium can 
be studied experimentally through systematic measure- 
ments of the pt spectra of open heavy flavor, i. e. hadrons 
composed of a heavy and a light quark. While the full re- 
construction of D meson decays at the Relativistic Heavy 
Ion Collider (RHIC) is reported for d+Au collisions [8|, 
indirect measurements of open heavy flavor via semi- 
leptonic decays are available for p+p and d+Au collisions 
at ^/J^T = 200 GeV 0,0, [3 as well as for Au+Au colli- 
sions at 130 and 200 GeV [nlE3 In P+P collisions, the 
extracted electron px spectrum from heavy flavor decays 
is in reasonable agreement with perturbative quantum 
chromodynamics (pQCD) calculations in next-to-leading 
order. However, the data leave room for contributions 
from further production mechanisms in which the heavy 
quarks are not created in the initial hard parton scatter- 
ing, e.g. via jet fragmentation [|J. In d+Au collisions, 
no indications for strong cold nuclear matter effects were 
found 0, El • F° r Au+Au collisions of different central- 
ity, the total electron yield from heavy flavor decays was 
observed to scale with the nuclear overlap integral (Taa) 
as expected for point-like pQCD processes |l2j. How- 
ever, these electrons show an azimuthal anisotropy with 
respect to the reaction plane |l3l | , consistent with the no- 
tion of charm quark flow in Au+Au collisions. It has been 
pointed out that if the charm quarks flow along with the 
bulk of the medium, this is evidence for thermalization 
of charm. In this situation, the medium modifications of 
the charm spectrum should be substantial 0] . 

In this Letter, we report on the px spectra of non- 
photonic electrons, (e + + e~)/2, measured at mid- 
rapidity (\rj\ < 0.35) up to px = 5 GeV/c by the PHENIX 
experiment in Au+Au collisions at ^/s^n = 200 GeV. 
The photonic electron background was removed by a 
cocktail subtraction, in contrast to the converter subtrac- 
tion used in 0] , where a subset of the current data sam- 
ple was analyzed. The converter method is better suited 
for a determination of the total yield of heavy flavor elec- 
trons, while the cocktail subtraction used here provides 



a precision measurement of the spectral shape. The nu- 
clear modification is then determined by comparing the 
spectra to those in p+p collisions . 

The data used in this analysis were collected by the 
PHENIX detector jl^ during the 2001 run of RHIC. A 
coincidence of the beam-beam counters (BBC) and the 
zero degree calorimeters (ZDC) provided the minimum 
bias trigger (92.2^3q % of the Au+Au inelastic cross sec- 
tion). The centrality was determined by the correlation 
between the multiplicity measured by the BBC and the 
energy of spectator neutrons measured by the ZDC. After 
restricting the vertex range to \z\ < 20 cm to eliminate 
background originating from the central magnet, a data 
sample of 25 x 10 6 minimum bias events was analyzed. 

For the electron analysis, charged particle tracks were 
reconstructed with the drift chamber and the first layer 
of pad chambers of the PHENIX east-arm spectrome- 
ter < 0.35, A(f> = 7r/2), as discussed in detail else- 
where |l2( . Tracks were confirmed by matching hits in 
the electromagnetic calorimeter (EMC) within 2 a in 
position. Electron candidates had at least three associ- 
ated hits in the ring imaging Cerenkov detector (RICH). 
After an additional cut on the correlation between the 
momentum p and the energy E deposited in the EMC 
(— 2cr < (E — p)/p < 3a), the only background remaining 
in the electron sample was due to accidental coincidences 
between RICH hits and hadron tracks. This background 
was estimated (~15 % at low px in central collisions, de- 
creasing towards high px and for peripheral events) and 
subtracted statistically by an event-mixing method. 

The raw electron spectra were corrected as a function 
of pt for geometrical acceptance and reconstruction ef- 
ficiency [l2j. The multiplicity dependent efficiency loss 
was estimated by embedding simulated electrons into real 
events. This loss does not depend on pt and increases 
from 5 to 26 % from peripheral to central collisions. The 
Id systematic uncertainty of all corrections is 11.8 %, 
after correction for the effect of finite bin width in px- 
The fully corrected inclusive electron spectrum is shown 
in Fig.^a) for minimum bias collisions. 

The spectra of electrons from heavy flavor decays were 
determined by subtracting cocktails of background con- 
tributions from other sources from the inclusive data. 
The most important background is the ir° Dalitz de- 
cay which was calculated individually for each centrality 
class with a hadron decay generator using parameteriza- 
tions of measured ir° Q and ^ 0] spectra as input. 
The spectral shapes of other light hadrons h were ob- 
tained from the pion spectra, assuming a universal spec- 
trum in rriT = -\/p T + n%i. Within this approach the 
ratios h/ir° are constant at high px with the values [llj| : 
77/71-° = 0.45 ± 0.10, p/ir° = 1.0 ± 0.3, uj/tt° = 1.0 ± 0.3, 
rj'/ir = 0.25 ± 0.08, and 0/tt = 0.40 ± 0.12. Only the 77 
contribution is of any practical relevance, and the chosen 
parameterization is in good agreement with the measured 
77 meson spectra |l8j . Another major electron source is 
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FIG. 1: (Color online) Inclusive and non-photonic elec- 
tron invariant yields in minimum bias Au+Au collisions at 
^/sjvjv = 200 GeV, compared with contributions from all 
background electron sources included in the cocktail (a). In- 
variant yields of electrons from heavy flavor decays for differ- 
ent Au+Au centrality classes, scaled by powers of ten for clar- 
ity. Curves are the best fit to the p+p reference scaled with 
the appropriate nuclear overlap integrals (Taa) (b). The error 
bars (brackets) correspond to statistical (systematic) uncer- 
tainties in both panels. 



the conversion of photons, mainly from ir° — > 77 decays, 
in material in the acceptance (ssl% X/Xq). The spectra 
of electrons from conversions and Dalitz decays are very 
similar. In a GEANT simulation of 7r° decays, the ratio of 
conversion electrons to Dalitz electrons was determined 
to be 1.25 ± 0.10, essentially px independent. Contribu- 
tions from photon conversions from other sources were 
taken into account as well. Electrons from kaon decays 
(K e ^), determined in a GEANT simulation based on mea- 
sured kaon spectra and electrons from external as 
well as internal conversions of direct photons [rHEoj l were 
included. 

All background sources are compared with the in- 
clusive data in Fig. QJa). Further background from 
J fil) — > e + e~ decays and from Drell-Yan pairs [2l| is 
negligible. A possible low mass dilepton enhancement 



through 7T + 7T — > p — > e + e~, as reported in Pb+Pb 
collisions at the SPS [22j, would constitute another back- 
ground source which is neglected here since the estimated 
p contribution in the absence of enhancement is small 
(< 1% at all pr). The total cocktail systematic uncer- 
tainty increases from 10 % (at px — 0.4 GeV/c ) to 
15 % (at px =5 GeV/c), dominated by the system- 
atic error of the pion input spectra (« 8-10 %). Other 
systematic uncertainties, mainly the tj/tt normalization 
and, at high px, the contribution from direct radiation, 
are much smaller. The background cocktail calculated 
here and the photonic electron background measured via 
the converter method 01 agree within 10 %. 

After subtracting the cocktail from the inclusive elec- 
tron data, the invariant spectrum of electrons from heavy 
flavor decays is shown in Fig.^a) for minimum bias colli- 
sions. For pt > 2 GeV/ c the signal to background ratio is 
larger than one. Fig.^b) shows the electron spectra from 
heavy flavor decays in four centrality classes, 0-10 %, 10- 
20 %, 20-40 %, and 40-60 % central collisions. More 
peripheral collisions have insufficient electron statistics 
to reach px = 5 GeV/c. 

PHENIX has also measured electrons from heavy flavor 
decays in p+p collisions at y/s = 200 GeV |IJ . The curves 
shown in Fig.^b) depict the best fit of the corresponding 
spectrum from p+p collisions, scaled by the nuclear over- 
lap integral (Taa) calculated within a Glauber model Q 
for each Au+Au centrality class. At low px the Au+Au 
spectra are in reasonable agreement with the p+p fit in 
all centrality bins, but a clear suppression of the spectra 
in Au+Au with respect to p+p develops towards high 
Pt- 

To quantify this effect we calculate for each individual 
bin in px the nuclear modification factor Raa defined as 



Raa = 



dN A 



u-\-Au 



(T. 



AA, 



x da 



(1) 



p+p 



where dNAu+Au is the differential electron yield from 
heavy flavor decays in Au+Au collisions and da p + p is 
the corresponding differential cross section in p+p colli- 
sions in any given px bin. 

Fig. |21 shows Raa as a function of px in the four 
Au+Au centrality classes. At lowpr, the electron Raa is 
consistent with one within substantial uncertainties in all 
centrality classes, in agreement with the observation of 
binary collision scaling of the total charm yield in Au+Au 
collisions at RHIC [l2| . Since the ratio of electrons from 
heavy flavor decays to background increases with increas- 
ing px, the systematic uncertainties of Raa decrease to- 
wards high px- Raa falls well below one for electron 
Pt > 2 GeV/c, providing clear evidence for heavy quark 
medium modifications. The observed high px suppres- 
sion is most significant for central collisions. However, 
the limited statistics do not allow to quantify the cen- 
trality dependence of heavy quark medium modifications. 
At the highest px, the electron Raa becomes as small as 
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FIG. 2: Nuclear modification factor Raa for electrons from 
heavy flavor decays as function of pr in Au+Au collisions at 
-y/ijvjv = 200 GeV for the different centrality classes. The 
error bars are statistical only. Error brackets (boxes) indi- 
cate the systematic errors related to the uncertainties in the 
Au+Au (p+p) measurements. The bands around one show 
the relative systematic uncertainties in Taa- For the most 
central collisions the 7r° Raa is shown for comparison pj . For 
these data, a 13 % pr independent systematic uncertainty 
(not plotted) represents the uncertainty in (Taa) and in the 
7T° yield normalization. 



that for 7T° 2], indicating substantial energy loss of heavy 
quarks in the medium. It is important to note that elec- 
trons at a given px originate from decays of higher pt D 
or B mesons, making model independent comparisons of 
Raa for light and heavy quarks impossible. 

The observed Raa is remarkable, as electrons with 
Pt > 3.5 GeV/c are expected to include significant con- 
tributions from B meson decays, and B mesons should 
suffer less than D mesons from medium modifications. 
Depending on their time scales, mechanisms by which 
heavy quarks are produced after the initial hard parton 
scattering, such as gluon splitting in jets, might lead to 
an attenuation at high px which then is due to a mix- 
ture of light parton and heavy quark energy loss in the 
medium created at RHIC. 

Fig-Blconfronts current model calculations [23H24] ] uti- 
lizing induced gluon radiation as the heavy quark en- 
ergy loss mechanism with the data for the 10 % most 
central collisions. The three curves (la-c) include elec- 
trons from charm decays only j23j. They correspond 
to different values of the time-averaged transport co- 
efficient q, which denotes the average squared trans- 
verse momentum transferred from a hard parton per 
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FIG. 3: Nuclear modification factor Raa for electrons from 
heavy quark decays as function of pr for the 10 % most central 
Au+Au collisions at ^snn = 200 GeV in comparison with 
predictions from models incorporating charm quark energy 
loss. The curves (la-c) and (2a-b) are taken from (23| and 
[24|. respectively, where contributions from B meson decays 
are included in (2a-b) only. Experimental uncertainties are 
shown as described in Fig. [5] 



unit path length while traversing the medium and, as 
such, is proportional to the density of scattering centers 
in the medium. Curve (la) applies for the case with- 
out the presence of any medium causing heavy quark 
energy loss (q — GeV 2 /fm). The q values of 4 and 
14 GeV 2 /fm, which correspond to the curves (lb) and 
(lc), lead to light quark energy losses which bracket 
the observed high pt suppression of neutral pions and 
charged hadrons. Predictions for charm energy loss 
from |23| for medium densities at the extreme high end 
of those allowed by the observed light quark energy loss 
are consistent with the electron data. Contributions from 
bottom decays, which are expected to be significant for 
Pt > 3 GeV/c, should lead to an increase of the pre- 
dicted Raa since b quarks arepresumably less affected 
by energy loss than c quarks 6]. The curves (2a-b) are 
taken from [24| ■ They include electrons from both D and 
B meson decays and correspond to initial gluon densities 
of dN g /dy = 1000(3500) for curve (2a(b)), respectively, 
which again lead to light parton energy losses bracketing 
the observed high pt pion suppression. However, at high 
Pt the predicted Raa for electrons from heavy flavor de- 
cays is larger than observed. The present data pose a 
challenge to existing calculations of radiative energy loss 
in the medium produced at RHIC, and will help to dis- 
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tinguish between different energy loss scenarios. 

In conclusion, we have measured electron spec- 
tra from heavy flavor decays in Au+Au collisions at 
y> r SNN = 200 GeV. In central collisions, nuclear modi- 
fication factors Raa << 1 are observed at high pr, pro- 
viding clear evidence for strong medium effects. Current 
models involving energy loss via induced gluon radiation 
for heavy quarks traversing the medium created in heavy 
ion collisions at RHIC are challenged by the data even 
considering extremely high medium densities. 

We thank the staff of the Collider-Accelerator and 
Physics Departments at BNL for their vital contribu- 
tions. We acknowledge support from the Department of 
Energy and NSF (U.S.A.), MEXT and JSPS (Japan), 
CNPq and FAPESP (Brazil), NSFC (China), CNRS- 
IN2P3 and CEA (France), BMBF, DA AD, and AvH 
(Germany), OTKA (Hungary), DAE and DST (India), 
ISF (Israel), KRF and CHEP (Korea), RMIST, RAS, and 
RMAE (Russia), VR and KAW (Sweden), U.S. CRDF 
for the FSU, US-Hungarian NSF-OTKA-MTA, and US- 
Israel BSF. 



* Deceased 

* PHENIX Spokesperson zajc@nevis.columbia.edu 

[1] K. Adcox et al, Phys. Rev. Lett. 88, 022301 (2002). 
[2] S.S. Adler et al, Phys. Rev. Lett. 91, 072301 (2003). 



[3] I. Arsene et al., Phys. Rev. Lett. 91, 072305 (2003). 
[4] B.B. Back et al, Phys. Lett. B578, 297 (2004). 
[5] J. Adams et al., Phys. Rev. Lett. 91, 172302 (2003). 
[6] Y.L. Dokshitzer and D.E. Kharzeev, Phys. Lett. B519, 
199 (2001). 

[7] J.-Y. Ollitrault, Phys. Rev. D46, 229 (1992); S. Voloshin 
and Y. Zhang, Z. Phys. C70, 665 (1996). 

[8] J. Adams et al., Phys. Rev. Lett. 94, 062301 (2005). 

[9] S.S. Adler et al, hep-ex/0508034 
[10] S. Kelly et al, J. Phys G30, S1189 (2004). 
[11] K. Adcox et al., Phys. Rev. Lett. 88, 192303 (2002). 
[12] S.S. Adler et al, Phys. Rev. Lett. 94, 082301 (2005). 
[13] S.S. Adler et al., Phys. Rev. C72, 024901 (2005). 
[14] G.D. Moore and D. Teaney, Phys. Rev. C71, 064904 
(2005). 

[15] K. Adcox et al, Nucl. Instrum. Methods A499, 469 
(2003). 

[16] GEANT 3.21, CERN program library. 

[17] S.S. Adler et al, Phys. Rev. C69, 034909 (2004). 

[18] S.S. Adler et al., in preparation. 

[19] L.E. Gordon and W. Vogelsang, Phys. Rev. D50, 1901 
(1994). 

[20] S.S. Adler et al, Phys. Rev. Lett. 94, 232301 (2005). 
[21] S. Gavin, P.L. McGaughey, P.V. Ruuskanen, and R. 

Vogt, Phys. Rev. C54, 2606 (1996). 
[22] G. Agakichiev et al, Phys. Lett. B422, 405 (1998). 
[23] N. Armesto, S. Dainese, C. Salgado, and U. Wiedemann, 

Phys. Rev. D71, 054027 (2005). 
[24] M. Djordjevic, M. Gyulassy, R. Vogt. and S. Wicks. 

|nucl-th/ 0507019 



